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Iodine is distinguished from other elements used in dark matter direct detection experiments both 
by its large mass as well as its large magnetic moment. Inelastic dark matter utilizes the large mass 
of iodine to allay tensions between the DAMA annual modulation signature and the null results 
from other experiments. We explore models of inelastic dark matter that also take advantage of 
the second distinct property of iodine, namely its large magnetic moment. In such models the cou- 
plings are augmented by magnetic, rather than merely electric, interactions. These models provide 
simple examples where the DAMA signal is compatible with all existing limits. We consider dipole 
moments for the WIMP, through conventional magnetism as well as "dark" magnetism, including 
both magnetic-magnetic and magnetic-electric scattering. We find XENONIOO and CRESST should 
generically see a signal, although suppressed compared with electric inelastic dark matter models, 
while KIMS should see a modulated signal comparable to or larger than that of DAMA. In a large 
portion of parameter space, de-excitation occurs promptly, producing a ~ 100 keV photon inside 
large xenon experiments alongside the nuclear recoil. This effect could be searched for, but if not 
properly considered may cause nuclear recoil events to fail standard cuts. 

PACS numbers: 12.60.Jv, 12.60. Cn, 12.60.Fr 



I. INTRODUCTION 

The DAMA annual modulation signature has persisted 
and grown to 8.9 a [T]. The constraints imposed by 
null results from many sources [2HS] would seem to ex- 
clude a dark matter interpretation of the signal within 
the framework of conventional weakly interacting mas- 
sive particles (WIMPs) scattering with spin-independent 
interactions. In this context, a wide range of propos- 
als have been put forward to explain the apparent ten- 
sion, such as spin-coupled WIMPs [Tllll). mirror mat- 
ter [121 [13] , resonant dark matter [2] , exothermic dark 
matter [EllTB], momentum-dependent scattering [TT/ISI, 
and light WIMPs with or without ion channeling 19 26J. 

Inelastic dark matter (iDM) [57] was originally pro- 
posed to explain DAMA not through a change in cou- 
plings, but rather through a change in kinematics. It 
brought into focus the large mass difference between io- 
dine (DAMA) and germanium (CDMS). The iDM frame- 
work has three basic elements 

• The existence of an excited WIMP state x* i with 
6 = 771^. — ra^ ^ ^^^v^ , where ^^x i^ the reduced 
WIMP-nucleus mass and v is the WIMP velocity. 

• An allowed inelastic scattering against the nucleus 

X + N — ^ X* + ^• 

• A forbidden, or highly suppressed elastic scattering 

X + N ^X + N. 

Inelastic scattering alters the minimum velocity required 
for a dark matter particle to scatter a nuclei N with recoil 
energy Er, 
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These kinematics change the spectrum of the signal 
(shifting it to higher recoil energy), favor heavy targets 
(such as DAMA's iodine) over lighter targets (such as 
CDMS's germanium), and enhance the modulated signal 
compared to the unmodulated part (favoring modulation 
experiments). Together, these features allow a signal at 
DAMA, while highly suppressing or eliminating the sig- 
nals at other experiments 28, 29 . 



In most models of iDM [271 150^^5] scattering arises 
dominantly through a vector current, such as a the Z- 
boson, or a new, light vector. However, many interac- 
tions can mediate inelastic transitions. For instance, re- 
cently 44J argued that a coupling to the proton's spin 
would significantly weaken constraints from existing ex- 
periments. While standard vector interactions are still 
viable, the strong constraints from ZEPLIN-III [33] and 
CRESST [5], in particular, require high rates of mod- 
ulation and departures from a Maxwellian halo [461 - 148] 
(although such departures are expected when sampling 
the tail of the distribution [iSllig] ). 



In this short letter, we point out that DAMA's iodine 
target also distinguishes itself from other detector's tar- 
gets through its magnetic properties. Iodine's large nu- 
clear magnetic moment compels us to consider models 
with dominantly magnetic interactions as a way to rec- 
oncile DAMA's positive detection claims with the other 
null results. We concentrate on dipole-dipole scattering, 
which can arise from conventional magnetism or a dipole 
of a dark force. In both cases, we shall see parameter 
space opens, with inelastic dipole-dipole scattering es- 
sentially unconstrained. 
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FIG. 1: The weighted-atomic mass and weighted-magnetic 
dipole moment (Eq. H in units of the nuclear magneton /ijv 
of various dark matter search targets. (C,0 and Ca,Ar have 
been shifted slightly so as not to overlay each other.) 



II. MAGNETIC INELASTIC DARK MATTER 

If one wants to understand how DAMA could have a 
positive signal while other experiments do not, there are 
many directions one can pursue. Narrowing the focus 
on nuclear recoils induced by WIMP coUisions, we must 
examine what the differences are between Nal and the 
other existing targets. 

The original iDM proposal focused on a single dimen- 
sion, namely the kinematical properties of iodine. As it 
is much heavier than many targets, in particular germa- 
nium, this allowed a significant departure from conven- 
tional elastic expectations. The fact that DAMA focuses 
on relatively high energies (~ 20-1- keV^ off iodine as- 
suming the standard quenching factor qj = 0.08) and 
modulation gave additional changes when comparing to 
elastic scattering limits, but ultimately the key distinc- 
tion was the kinematical change of a heavy target. 

This simple one-dimensional analysis is important, but 
iodine's magnetic properties also distinguish it from most 
other target nuclei. The quantity that we will see is most 
relevant is the weighted dipole moment 
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where fi, ^i, and Si are the elemental abundance, nu- 
clear magnetic moment, and spin, respectively, of isotope 
i. We show in Fig. [T] the abundance- weighted atomic 
masses, and the weighted dipole moment of various tar- 
get nuclei. We see that while tungsten (W) has a large 
mass, its magnetic moment is rather small. Fluorine (F) 
and sodium (Na) have large magnetic dipoles but are very 
light. Xenon (Xe) has a couple of isotopes with apprecia- 
ble dipoles, however, they are insufficient to make it com- 
petitive with iodine. The combination of large mass and 



large dipole makes the iodine target used by DAMA quite 
unique among the nuclear targets, with only KIMS' j50] 
cesium (Cs) target similar in its qualitative features. The 
iodine dipole arises dominantly from the angular momen- 
tum of unpaired protons |51) , with additional contribu- 
tions from the neutron and proton spin. 

We are therefore led to consider models that make both 
kinematical and magnetic distinctions between targets. 
Since its proposal, the focus of iDM model building has 
dominantly been on electrically coupled WIMPs, either 
directly to charge, or to some combination of the mass 
number A and the atomic number Z, such as through the 
Z°-boson. Since we wish to take advantage of the large 
magnetic dipole of iodine, we instead focus on models of 
magnetically-coupled inelastic dark matter (MiDA-I). 



III. SCENARIOS FOR MIDM 

The magnetic interactions of a WIMP can appear at 
different orders in the multipole expansion. The first 
order, namely a magnetic monopole, is interesting but 
problematic ^\. Instead we choose to focus on the case 
of a magnetic dipole which has a sizable interaction with 
the magnetic dipole of iodine. However, a magnetically 
interacting WIMP also feels a velocity-suppressed inter- 
action with the charge of the nucleus, thus one cannot 
simply consider scattering off magnetic moments. For 
iodine the contribution from Z^v^ is subdominant to 
fj?, but for magnetically-challenged nuclei, such as W, 
or even Xe, the charge coupling can dominate the scat- 
tering. 



A. Dipole-Dipole Inelastic Scattering 

The idea that the WIMP could have a magnetic dipole 
has been long studied (see., e.g., [55H55] .) The dipole 
operator is naturally off-diagonal HH [SI] , and mediates 
transitions between the ground state x ^nd the excited 
state X* 1 
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where /i^ is the dipole strength and tr^^ — i[jf^,j^]/2. 
[55] considered such transitions in the early universe for 
dark matter in the range of few keV — few MeV. [BU] con- 
sidered inelastic WIMP dipole-nuclear charge scattering 
to explain DAMA. Such an interaction, however, does 
not significantly change the relative strength of the vari- 
ous experiments compared with charge-charge (i.e., vec- 
tor current) interactions, and the viability of the scenario 
found in |60| was largely because the significant con- 
straints from the CRESST experiment were ignored. [33] 
considered a related idea, studying the parameter space 
under the assumption of an iDM that couples to pro- 
ton nuclear spin exclusively, although no particle physics 
model generating the required interaction was found. 



Here we study inelastic dipole-dipole (DD) interac- 
tions, which occur in addition to dipole-charge interac- 
tions, and we shall see that this opens up significant pa- 
rameter space. In fact, for WIMP-iodine scattering DD 
cross section dominates over the dipole-charge (DZ) scat- 
tering. DZ scattering, on the other hand, may dominate 
the scattering in other targets (such as tungsten) and pro- 
vide signals at other experiments. The direct detection 
scattering rate is given by. 
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where F'^[Er] is the usual nuclear form- factor and 
F'^[Er] is the nuclear magnetic dipole form-factor. In 
the above expressions we have kept the leading terms in 
the expansion of Er, S, and v. We note the inelastic en- 
hancement of the destructive interference in the DZ cross 
section which leads to a further reduction of the signal in 
experiments which utilize a magnetically-challenged tar- 
gets. We have checked that our formulae agree in the 
elastic limit with those in [521 EI] ■ 

We show the allowed parameter space in Fig. [2] To 
determine these parameter ranges, we use the binned 
DAMA data of [1] , taking the 2-8 keVee (upper row) and 
2-14 keVee (lower row) bins and calculating a x^ param- 
eter for each point in {m-^, S, fj,^) parameter space, using 
a quenching factor qj = 0.08. We have neglected scatter- 
ing off Na, as this is only important for 5 < 30 keV. For 
the dark matter velocity distribution, we take a Maxwell- 
Boltzmann distribution with parameters vq = 220 km/s 
and Vesc = 550 km/s. We find the global x^ mini- 
mum and plot slices of allowed parameter space at con- 
stant m,, showing the 90% and 99% confidence regions 
(Ax^ — 6.25 and 11.34, respectively), after ensuring 
that the best fit point gives a good fit to the data (i.e., 
x'^/dof < 1). We note that higher mass dark matter 
are disfavored in fits that include the higher energy bins. 
The high recoil energy region strongly depends on the 
behavior of the form-factor and the two rows therefore 
serve as an illustration of the effects of form-factor un- 
certainties. We discuss the nuclear form factor uncertain- 
ties further below. We see that the scenario of inelastic 
dipole-dipole scattering is essentially unconstrained for a 
wide range of the inelasticity 5. The allowed range of 
H^ ^ few X 10~"^/iN is comparable to what could arise 



from an electroweak loop, for instance, or if the WIMP 
is composed of charged constituents, bound by a force 
confining at the TeV scale, such as a technibaryon ^ . 

For our constraints, we use the full CDMS dataset. 
For ZEPLIN-III [H2] and XENONIO ^ we use re- 
cently published reanalyses focused on inelastic dark 
matter. The experiments closest to constraining the 
allowed DAMA region in our scenario are KIMS and 
ZEPLIN-III; ZEPLIN-III and CRESST have weakened 
limits relative to standard iDM 46J, due to their sup- 
pressed magnetic couplings. For Zeplin-III, a 90% CL 
limit of 5.4 signal events |S2] was extracted using the 
maximum patch technique. However, it is worth noting 
that the maximum patch technique used in [62j would 
likely have different sensitivity for MiDM's altered re- 
coil spectrum. Finally, we note that other lighter targets 
with large dipole moments, such as flourine and sodium, 
play essentially no role in the scenarios we consider. This 
is so because the inelastic threshold, S, is too large com- 
pared with the kinetic energy in the center of mass frame, 
/^A'x^^/2, even for the highest speeds available in the 
halo. 

Following the scattering against the nucleus, the ex- 
cited DM state leaves the target and de-excites shortly 
after with a lifetime given by r~^ = ii^S^/ir. For a 
typical splitting of (5 = 120 keV and a dipole moment 
/i, = 3 X 10~'^^Ar the lifetime is t = 5.1 ^s, leading 
to an average decay length of ur ~ 1.5 m. This can 
give de-excitations inside the detector for a fraction of 
the events depositing a photon of energy S ^ 100 keV. 
We note that in more complicated scenarios with mul- 
tiple states, the de-excitation may occur to a different 
state other than the ground state, so the photon energy 
needn't be precisely equal to that of the excitation en- 
ergy. (A similar phenomenology can occur in resonant 
dark matter (rDM) , although with a photon of typically 
much higher energy, ~ 500 keV -1 MeV [14].) 

Thus, the decay length of the excited state can be 
within an order of magnitude of the size of modern de- 
tectors. Consequently, additional light could be added to 
WIMP scatters, making them appear either anomalous 
or more similar to electron recoils. In particular, in a 
large xenon detector (such as XENONIOO or LUX) such 
an effect could potentially push the WIMP signal up, 
out of the nuclear recoil band, while leaving a signal at 
small detectors. Such decays could either give a signal of 
this scenario (if a second population appears, in addition 
to the primary one) , or remove most of the nuclear recoil 
band signal, if the majority of the decays occur inside the 
detector. At DAMA in particular, the majority of the de- 
cays would occur either outside the experiment or inside 
a different crystal from the one where the initial scatter 
occured. Events separated by > 600 ns are not counted 
in the multi-hit rate, and consequently existing cuts will 
not veto against these events. This higher energy modu- 
lation signal at S, due to the de-excitation, is consistent 
with their observation of no modulation above 90 keVee 
[I], but could be searched for with finer binning. Other 
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FIG. 2: Allowed ranges of parameter space for inelastic scattering of a magnetic dipole (dominated by dipole-dipole scattering 
for iodine). Purple (dark) and blue (light) regions are 90% and 99% allowed confidence intervals, respectively, for a x^ fit to the 
DAMA modulation. Constraints are: CDMS (solid), CRESST-II (short light dashed), XENONIO (short dark dashed), KIMS 
(long dark dashed), and ZEPLIN-III (long light dashed). In the upper row, we utilized only the 2 — 8 keVee bins in DAMA, 
and in the lower row we used the entire range of 2 — 14 keVee. This serves to illustrate the strong dependence of the allowed 
DAMA parameter space for heavy dark matter on the dipole form factor's behavior at high Er as discussed in Sec. [Ill C] in 
detail. 



experiments, however, and in particular, xenon experi- 
ments with large volumes might reject these events by 
vetoes designed to remove multiple scatter events since 
the de-excitations occur displaced with respect to the 
initial scatter and yield a distorted signal shape. These 
de-excitations are a potentially exciting and useful signa- 
ture, and we defer additional discussion to future work. 



B. Dark Dipole Scattering 

While a sizeable dipole of standard magnetism is possi- 
ble, the WIMP could instead have a large dark magnetic 
dipole under some additional U{1) which is kinetically 
mixed with electromagnetism through eFfj^i^F'^ /2 [51] . 
Such a feature would be particularly natural in compos- 
ite inelastic models (CiDM) [351 |5S]. While the interac- 
tion is suppressed by the mixing parameter e, it can be 
enhanced in that the constraints on new dark charged 
constituents are far weaker than those for ordinary elec- 
tromagnetism. 

In the limit that the mediator is very light m\, <C 9^ = 
2mffEii, the parameter space in Fig.[2]is unchanged, but 
with fiy^ — >■ /i^e. We can also consider m^, ^ q^. In this 



case, the cross section is a simple modification of the 
QED formula, Eq. ^, 
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For the case where ttia' = 70 MeV, the parameter space 
is shown in Fig.[3j In this scenario, neglecting Na scatter- 
ing at DAMA is an excellent approximation, as sodium 
is additionally suppressed relative to iodine by the g^ 
factors. An interesting effect from the additional g^ in 
the numerator is that the DAMA fit improves for smaller 
values of S, with a closer fit to the modulation shape. 
Larger values of the mediator mass, uia' , result in an in- 
crease (decrease) in the allowed parameter space shown 
in Fig. Is] for small (large) WIMP masses, rriy^ < 70 GeV 
("'•x ^ 70 GeV). The cross sections required to fit 
DAMA increase for a heavier mediator, with a shift to 
lower values of S. 



C. Form Factor Uncertainties 

As is well-known, spin-independent elastic scattering is 
suppressed at higher energies by the nuclear form factor. 
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FIG. 3: Allowed ranges of parameter space for dark dipole scattering -with mediator mass jn^/ = 70 MeV. In this plot, we 
have set the dark gauge coupling ad = a. Lines are labeled as in Fig. p] In the upper row, we utilized only the 2 — 8 keVee 
bins in DAMA, and in the lower row we used the entire range of 2 — 14 keVee. This serves to illustrate the strong dependence 
of the allowed DAMA parameter space for heavy dark matter on the dipole form factor's behavior at high Er as discussed in 
Sec. HITCl in detail. 



Explicit measurements of F'^[Eb\ have been made (see 
the discussion in |66j). and are in good agreement with 
the Helm form factor, which we employ for the contribu- 
tions from dipole-charge scattering. 

In contrast, there are no explicit measurements of the 
magnetic dipole form factor, or well-developed calcula- 
tions available in the literature. The magnetic dipole 
moment arises from a variety of contributions in different 
nuclei, with a significant contribution often from proton 
angular momentum (in iodine, in particular) as well as 
from proton and neutron spin [51 . Odd group and sin- 
gle particle models give soft form factors for the magnetic 
moment [ST]. [53] employs an estimated "spin radius", 
which gives a similarly soft form factor. However, these 
approaches neglect the nuclear spin contribution, whose 
form factors tend to have a hard tail, and the contribu- 
tion of this should be taken into account. 

Thus, we employ an approximation for the magnetic 
form factor: we take the contributions of angular mo- 
mentum for iodine and cesium to be arising from, re- 
spectively, a 2d-shell orbital and a Ig-shell orbital, as ex- 
pected in the nuclear shell model (see for e.g. [17]), which 
we calculate explicitly using a SHO model of the nucleus. 
For the SHO, we use a harmonic oscillator parameter of 
h = 2.282 fm, which characterizes the size for an iodine 
nucleus ^51^. We add to this a contribution from nuclear 



spin, which is taken from [5T] (I) and [SS] (Cs). Since 
only squared-form factors are available, we take a square 
root and assume no sign change takes place in order to 
combine with the angular momentum contribution. We 
weight the angular momentum and spin contributions to 
the magnetic dipole moment as suggested by the values 
in [51] [68] • Specifically, we use 
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which characterize the contributions of the proton spin 
and angular momentum to the magnetic dipole moment 
for iodine and cesium which are both proton-odd. The 
predicted scattering rate will depend on these coefficients 
in the combination, but of these the largest dependence 
is in cesium due to its destructive interference. We have 
chosen the more conservative result that comes from 
using renormalized g factors |68) . It is worth noting 
that only after using these renormalized factors, is the 
magnetic dipole moment of Cs accurately predicted in 
[68) . Using the unrenormalized factors would lead to a 



stronger KIMS limit, especially at higher WIMP masses. 
The other heavy targets we consider have dipoles dom- 
inated by neutron spin, so we use the standard neutron 
spin form factor for their scattering. 

We should emphasize that combining the spin and an- 
gular momentum contributions in this way is conserva- 
tive. Softer form factors tend to give better agreement 
for higher masses, and the destructive interference we 
assume for Cs means that the overall magnetic dipole 
form factor falls off slowly with En, since the cancella- 
tion ceases to be effective as the angular momentum con- 
tribution drops off rapidly. This results in particularly 
stringent limits from KIMS. Given these uncertainties, 
we have employed an overall event rate constraint from 
KIMS, requiring the event rate from 3 — 8 keVee (equiva- 
lently 20 — 100 keVR,) to be less than the 90% confidence 
level region of the observed rate .28 ± .16 cpd/kg. Note 
that this energy recoil range roughly corresponds to the 
2 — 8 keVee region we analyze for DAMA. On the other 
hand, since signals from Xe and W are dominated by 
dipole-charge scattering, their limits should be relatively 
robust. 

We emphasize that our choice of a DAMA x^ fit to 
the 2-8 keVee bins shown in the first row of Figs. [2] and 
[3] is insensitive to the form factor at Eft > 100 kcVR. 
This row effectively illustrates the parameter space in 
case of a softer form-factor, since the 8 — 14 keVee bins 
would be less important for such a form factor. If the 
higher energy bins are included, as shown in the lower 
row of each figure, then a harder form-factor tail, such 
as the spin form factor, results in a poorer fit to DAMA 
for heavy dark matter masses (and in particular at large 
i5). For example, we can see in Fig. [2] and l3J how the 
DAMA region shrinks entirely when the higher energy 
bins are included for the m^ = 300 GeV plots, while 
the plots for rn^ = 70 GeV are nearly unaffected. Thus, 
whether or not heavy WIMPs are allowed in this scenario 
is very sensitive to the high energy behavior of the dipole 
form factor. Consequently, numerical calculations of the 
form factor would help to pin down this parameter space. 
Note that in the cases where the WIMP is composite, 
additional form factors are possible which may soften the 
spectrum at high energy, adding additional uncertainty 
beyond those of the nuclear magnetic dipole form factor. 



the DAMA modulation data is possible without violating 
constraints imposed by other searches. Given this largely 
untested scenario, we strongly encourage calculations of 
the nuclear magnetic dipole form factor, as this remain- 
ing uncertainty is important to determine the parameter 
ranges of interest. 

Magnetic-like interactions also change our expecta- 
tions for upcoming experiments. When the scattering at 
DAMA is enhanced by a sizeable dipole-dipole contribu- 
tion, the signal experiments with small magnetic dipoles 
is suppressed by a comparable amount. At any of the 
xenon experiments, CRESST or CDMS, this is an ad- 
ditional suppression of a factor of a few. Moreover, as 
the lifetime of the excited state can be short, decays of 
the excited state in the detector can mask it as a dou- 
ble scattering event, potentially suppressing signal even 
farther in large xenon experiments. Thus, compared to 
standard iDM, a smaller signal would be expected in the 
current run at XENONIOO. On the other hand, KIMS 
with its Csl target retains sensitivity to this scenario 
and could search for it by either repeating their anal- 
ysis with additional exposure or looking for modulation. 
Of course, should the WIMP have both electric and mag- 
netic dipoles mediating the inelastic transitions, then in- 
termediate scenarios are possible. However, in all cases 
it seems that with ~ 4000 kg day of summer exposure, 
XENONIOO should be able to definitively test any sce- 
nario in which scattering of the magnetic dipole plays 
a major role. Because de-excitation can occur within 
the detector volume, standard cuts may miss a sizable 
fraction of these events. However, if attention is paid 
to double-scatter coincident events, xenon experiments 
should be capable of detecting both the inelastic WIMP 
signal and its magnetic origin. 

In conclusion, we have considered models of iDM where 
dark matter has a (EM or dark) magnetic dipole, which 
plays a major role in its scattering off of nuclei. The 
presence of the large magnetic dipole moment of iodine 
in addition to its large mass widens the parameter space 
where the modulation signal at DAMA is consistent with 
the null results from all other experiments. All mod- 
els predict some signal at the XENONIOO, KIMS, and 
CRESST experiments, with the usual iDM spectral fea- 
tures, and a large modulation amplitude. With new data 
on the horizon, these models should be tested soon. 



IV. CONCLUSIONS 

The iodine target used at DAMA is set apart from 
other direct detection targets by two intrinsic properties: 
its large mass and nuclear magnetic moment. Therefore, 
inelastic magnetic transitions which naturally utilize sig- 
nificant dipole-dipole interactions in addition to dipole- 
charge can explain the discrepancy between DAMA and 
other searches. In this letter, we considered this pos- 
sibility in detail for WIMPs with magnetic dipoles or 
alternatively dark magnetic dipoles. We find that large 
parts of parameter space are available where a good fit to 



Acknowledgments 

We would like to thank Laura Baudis, Pierluigi Belli, 
Rafael Lang, Dan McKinsey, and Peter Sorensen for use- 
ful discussions pertaining to the different experiments. 
We are grateful to Wick Haxton and David Dean for 
their insightful and extremely patient discussions of the 
spin structure of nucleons. We thank T.Banks and S. 
Thomas for useful discussions. They have informed us 
that they are writing a paper with J-F. Fortin on the 
signals for dark matter with electromagnetic dipoles in 



terrestrial detectors. While there is some overlap with 
our work (they also emphasize the importance of the 
large magnetic moment of Iodine), their models do not 
have inelastic scattering and differ significantly from our 



own. SC is supported under DOE Grant #DE-FG02- 
91ER40674. NW is supported by DOE OJI grant #DE- 
FG02-06ER41417 and NSF grant #0947827. lY is sup- 
ported by the James Arthur fellowship. 



[1] R. Bernabei et al. 
S10052-010-1303-9 



J2 



I3 

^. 

[6: 

[7 

[8 

[9 

[10 

[11 
[12 
[13 

[14; 

[15 

[16: 

[17 

[is: 

[19 
[20 

[21 
[22: 
[23 
[24 
[25 
[26 
[27 



(2010), doi:10.1140/epjc/ 

arXiv:1002.1028 [astro-ph.GA]^ 



E. Aprile et al. 
[astro-ph.COi] 



(XENONIOO), (2010), arXiv:1005.0380| 



Z. Ahmed et al. (The CDMS-II), (2009), |arXiv:0912. 3592 

[astro-ph.CO] . 

V. Sanglard (for the EDELWEISS), J. Phys. Conf. Ser., 

203, 012037 (2010), ,arXiv:0912.1196 [astro-ph.CO], . 

R. F. Lang an d W. Seidel, New J. P hys., 11, 105017 

(2009), arXiv:0906.3290 [astr o-ph.IM]| . 

V. N. Lebedenko et al. (ZEP LIN-IIlt; Phys. R ev. Lett., 

103, 151302 (2009), arXiv:0901.4348 [hep-ex] . 

P. Ulho M. Ka mionkowski, and P. Vogel, JHEP, 07, 044 

(2001), |arXiv:h ep-ph/0010036 . 

P. Belli, R. Cerulli, N. Fornengo, and S. Scopel, Phys. 

Rev., D66, 043503 (2002), arXiv:hep-ph/0203242 . 

C. Savage, P. Gondolo, and K. Freese, Phys. Rev., DTO, 

123513 (2004), arXiv:astro-ph/0408346 . 

R. Bernabei et al, Int. J. Mod. Phys., D13, 2127 (2004), 

arXiv :astro-ph/05014T2l . 

'MTTa irbairn and T. Schwetz, JCAP, 0901, 037 (2009), 

|arXiv:0808.0704 [hep-ph], . 

R. Foot, Mod. Phys. Lett., A19, 1841 



(2004), 



[28 
[29 
[30 

[31 



[35 
[36 

[37 
[38: 



D78, 043529 (2008), f 



|arXiv: astro-ph/0405362 

R^ Foot, Phys. Rev., 

|arXiv :0804.4518 [hep-ph] . [39 

Y. B ai and P. J. Fox, JHEP, 11, 052 (2009), 
|arXiv :0909.2900 [hep-ph] . [40 

R. Essig, J. Kaplan, P. Schuster, and N. Toro, (2010), 
arXiv:1004.0691 [hep-ph] . [41 

P. W. Graham, R. Harnik, S. Rajendran, and 
P. Saraswat, (2010), arXiv: 1004.0937 [hep-ph] . [42 

B. Feldstein, A. L. Fitzpatrick, and E. Katz, JCAP, 
1001, 020 (2010), a rXiv:0908 .2991 [hep-ph] . [43 
S. Chang, A. Pierce, and N. Weiner, JCAP, 1001, 006 
(2010), arXiv:0908.3192 [hep-ph] . [44 
P. BeUi et al, Phys. Rev., D61, 023512 (2000), 
arXiv:hep-ph/9903501 . [45 
A. Bottino, F. Donato, N. Fornengo, and S. Scopel, Phys. 

Rev., D69, 037302 (2004), |arXiv:hep-p h/0307303 . [46 

R. Be rnabei et al, Eur. Phys. J., C53, 205 (2008), 
[arXi^ :0710.0288 [astro-p h]]. [47 

A. Bottino, F. Donato, N. Fornengo, and S. Scopel, Phys. 
Rev., D78, 083520 (2008), ar Xiv:0806.4099 [hep -ph] . 
S. Chang, A. Pierce, and N. Weiner, (2008), 
arXiv:0808.0196 [hep-ph] . [49 

C. Savage, G. Gelmini, P. Gondolo, and K. Freese, 
JCAP, 0904, 010 (2009), arXiv:0808.3607 [astro-ph]'. [50 
F. Petriello and K. M. Zurek, JHEP, 09, 047 (2008), 
arXiv:0806.3989 [hep-ph] . [51 

D. Hooper, J. I. Collar, J. Hall, and D. McKinsey, 
(2010), ar Xiv: 1007. 10 05 [hep-ph] . [52 
D. Tucker-Smith and N. Weiner, Phys. Rev., D64, 



S. Chang, G. D. Kribs, D. Tucker-Smith, and N. Weiner, 

(2008), arXiv:080 7.2250 [hep-ph] _ . 

J. March-RusseU , C. McCabe, and M. McCuUough, 

(2008), arXiv:08 12.1931 [a stro-p"h]1 . 

N. Arkani-Hamed, D. P. Finkbeiner, T. R. Slatyer, 

and N. Weiner, Phys. Rev., D79, 015014 (2009), 

arXiv:0810.0713 [hep-ph] . 

M. Baumgart, C. Cheung, J. T. Ruderman, L.-T. Wang, 

and I. Y avin, JHEP, 04, 014 (2009), , arXiv:0901.0283 , 

[hep-ph] I . 



Y. Cui, D. E. Morrissey, D. Poland, and L. Randall, 
(2009), arXiv:0901.0557 [hep-ph] |. 

C. Cheung, J. T. Ruderman, L.-T. Wang, and I. Yavin, 
Phys. Rev., D80, 035008 (2009), _arXiv:0902. 3246 [hep-^ 

A. Katz and R. Sundrum, JHEP, 06, 003 (2009), 



arXiv:0 902j2ri_|hep-ph]| . 

D. P. Finkbeiner, T. R. Slatyer, N. Weiner, and I. Yavin, 

JCAP, 0909, 037 (2009), arXiv:0903.1037 [hep-ph] . 

D. S. M. Alves, S. R. Behbahani, P. Schuster, and J. G. 

Wacker, (2009), arXiv:090 3.3945 [hep-ph]. . 

D. E. Morrissey D. Pola nd, and K. M 'Zurek, JHEP, 

07, 050 (2009), arXiv:0904.2567 [hep-ph]i. 

M. M. Anber, U. Aydemir, J. F. Donoghue, and P. Pais, 

Phys. Rev., D80, 015012 (2009), |arXiv:0905.4260 [hep^ 

■phU- 

C. Arina, F.-S. Ling, and M. H. G. Tytgat, JCAP, 0910, 

018 (2009), arXiv:0907.0430 [hep-ph] . 

F. Chen, J. M. Cline, and A. R. Frey Phys. Rev., D80, 

083516 (2009), arXiv:0907. 4746 [hep-ph] | . 

D. E. Kaplan, G. Z. Krnjaic, K. R. Rehermann, and 
C.M.Wells, (2009), arXiv:0909.0753 [hep-ph],. 

A. Kum ar, D. Tucker-Sm ith, and N. Weiner, (2009), 



arXiv:0 910.2475 [hep-ph 

S 



and N. Weiner, (2010), 



Cha ng, R. F. Lang, 
| arXiv: 1007.2688 [hep-ph] . 

J.^Kopp, T. Schwetz, and J. Zupan, JCAP, 1002, 014 
(2010), arXiv:0912.4264 [hep-ph] . 

M. I. Lopes (ZEPLIN-IH), J. Phys. Conf. Ser., 203, 
012025 (2010). 

M. Kuhlen e t al, JCAP , 1002, 030 (2010), 

|arXiv:0912.2358 [astro-ph.GA] . 



R. F . Lang and N. Weiner, (2010), ar Xiv:1003.3664 [hep-| 
M. Alves, M. Lisanti, and J. G. Wacker, (2010), 



ph] 
D. S. 



|arXiv:1005.5421 [hep-ph 



ph 



M. Vogelsberger et al, (2008), |arXiv:0812.0362 [astro-| 
Phys. Conf. Ser., 120, 042021 



043502 (2001), arXiv:hep-ph/0101138 



STK. Kim (KIMS), J. 

(2008). 

M. T. ResseU and D. J. Dean, Phys. Rev., C56, 535 

(1997),,arXiv:hep-ph/9702290 . 

The flux of Standard Model monopoles is strongly 

bounded [? ]. A monopole under some other, broken 

(7(1) has to be confined, which complicates the analysis. 



and we leave it for future work. 
[53] J. Bagnasco, M. Dine, a nd S. D. Thomas, Phys. Lett., 

B320, 99 (1994), |arXiv: hep-ph/9310290 . 
[54] M. Pospelov and T. ter Veldhuis, Phys. Lett., B480, 181 

(2000), arXiv;hep-ph/0003010 . 
[55] K. Sigurdson, M. Doran, A. Kurylov, R. R. Caldwell, 

and M. Kamio nkowski, Phys. Rev., D70, 083501 (2004), 

|arXiv:astro-'ph/ 0406355 . 
[56] S. Gardner, Phys. R ev., DT9, 055007 (2009), 

|arXiv:0811.0967 [hep-ph]| . 
[57] W. S. Cho, J.-H. Huh, L-W. Kim, J. E. Kim, and 

B. Kyae , Phys. Lett., B687, 6 (2010), arXiv;1001.0579 
r [hep-ph[ . 
[58] H. An, S.-L. Chen, R. N. Mohapatra, S. Nussinov, and 

Y. Zhang, (2010), arXiv: 1004.3296 [hep-ph] . 
[59] H. K. Dreiner, H. E. Haber, and S. P. Martin, (2008), 

|arXiv :0812.1594 [hep-ph] . 
[60] E. Masso, S. Mohanty, and S. Rao, Phys. Rev., D80, 

036009 (2009) , arXiv:0906. 1979 [hep-ph] . 



[61] V. Barger, W.-Y. Keung, and D. Marfatia, (2010), 

arXiv: 1007.4345 [hep-ph] . ^ 

[62] D. Y. Akimov et al, (2010), |arXiv:1003. 5626 [hep-ex] 



[63] X. Collaboration, (2009), arXiv:0910.3698 [astro-ph.CO]] 
[64] B. Holdom, Phys. Lett., B166, 196 (1986). 



[65] M. Lisan ti and J. G. Wacker, (2009), arXiv:091 1.4483 
[hep-ph' 



[66] G. Dud a, A. Kem per, and P. Gondolo, JCAP, 0704, 012 

(2007), arXiv:he p-ph/0608035 . 
[67] J. R. EllTs^lmd R. A. Flores, Phys. Lett., B263, 259 

(1991). 
[68] P. Toivanen, M. Kortelainen, J. Suhonen, and J. Toiva- 

nen, Phys. Rev., C79, 044302 (2009). 
[69] The flux of Standard Model monopoles is strongly 

bounded [? ]. A monopole under some other, broken 

U{1) has to be confined, which complicates the analysis, 

and we leave it for future work. 



